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Abstract 

Nucleon-nucleon (NN) forces from chiral perturbation theory at next-to-leading (NLO) and next- 
to-next-to-leading order (NNLO) are applied to systems with two, three and four nucleons. At 
NNLO, we consider two versions of the chiral potential which differ in the strength of the two- 
pion-exchange (TPE) but describe two nucleon observables equally well. The NNLO potential 
leads to unphysical deeply bound states in the low partial waves and effects of the 3N forces, which 
appear first at this order, are expected to be large. We provide arguments for a reduction of the 
TPE potential and introduce the NNLO* version of the NN forces. We calculate nd scattering 
observables as well as various properties of and ''He with the NNLO* potential and find good 
agreement with the data and with predictions based upon the standard high-precision potentials. 
We find an improved description of the '^H and ^He binding energies. 
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1 Introduction 



Nuclear forces are derived in the chiral effective field theory approach in terms of an expansion in 
powers of where Q corresponds to a generic external momentum of nucleons and represents 

the typical hadronic scale (scale of chiral symmetry breaking) of the order of 1 GeV. That ratio is less 
than one if one considers processes with sufficiently low external momenta of the nucleons. In order 
to exclude contributions of high-momentum components in intermediate states, the nucleon-nucleon 
(NN) potential is multiplied by a regulator, which suppresses momenta larger than a certain cut-off 
A [1]. The latter has to be chosen below the scale A^7 The cut-off A should also not be taken too 
small in order not to suppress the relevant physics. The various coupling constants depend on the 
cut-off A in a way to compensate the changes in the low-energy observables induced by varying A. 
The remaining cut-off dependence of the observables can be removed by adding higher order terms 
to the effective potential [1]. Assuming naturalness for the various renormalized coupling constants 
in the underlying Lagrangian one can expect that contributions to the NN forces corresponding to 
higher powers v of the chiral expansion will decrease. This sort of nuclear interactions based on the 
most general chiral invariant effective Lagrangian formed out of pion and nucleon fields has been first 
proposed in [4] and formulated in detail in [5]. We followed a similar path, however extracting the 
nuclear forces from the Lagrangian in a different way. We refer to [6] where two and three- nucleon 
potentials have been derived using the method of unitary transformation. That method leads to 
energy independent and hermitean nuclear forces which are better suited for applications to systems 
with A > 2 than energy-dependent forces derived in old fashioned time-ordered perturbation theory 
like in [5]. In [7] we applied the forces at next-to-leading order (NLO), corresponding to (the counting 
index) u = 2, to the 3N and 4N systems. At this order NN phase shifts can be described only at 
rather low energies and only modestly. Nevertheless 3N and 4N binding energies were found to be 
within the same range as the ones found with high precision modern NN forces and also nd elastic 
and break-up observables at very low energies are similar to predictions generated by conventional 
forces. At that order the experimental nucleon analyzing power Ay is fairly well reproduced, which 
for conventional NN forces poses a serious puzzle [8]. This result, however, has to be considered as an 
intermediate step, corresponding just to NLO, where the ^Pj NN phase shifts could not be reproduced 
with sufficient accuracy. It is now of strong interest to explore the chiral forces in 3N and 4N systems 
at ncxt-to-ncxt-to-leading order (NNLO) corresponding to = 3 where the NN phase shifts are 
better reproduced. For the convenience of the reader we review briefly the NN forces in LO (i^ = 0), 
NLO and NNLO in section 2. 

It has been already pointed out in [9] that the strong central attraction caused by the numerically 
large values of the LEC's ci, C3, and C4 as determined in a analysis of vrN scattering leads to 
spurious deeply bound states in various two-nucleon angular momentum states. Though this has no 
observable consequences in the NN system within the realm of validity of the theory it is technically 
somewhat disturbing in treating 3N and 4N systems. Also ignoring 3N forces, which occur at NNLO 
the first time, and exploring only the NNLO NN forces leads to strong deviations from 3N data as we 
will show. It has to be expected that this will be remedied by including the NNLO 3N forces, which 
necessarily have to be taken into account at that order. Various consequences of the large values of 
the Ci's as well as the current situation in relation to the determination of the Cj's from other processes 
(such as ttN scattering) are discussed in section 3. Motivated by the findings of the boson-exchange 
(BE) models of the nucleon-nucleon interaction, we constructed the NNLO* potential by removing 
the A content from the LEC's C3 and C4 and refitting the contact interactions. The new values of 
the Cj's resulting from subtracting the A contributions lead to the NNLO* potential which is free of 

^In some cases it turns out to be possible to perform standard renormalization of the theory by taking the cut-off A 
to infinity [2, 3]. 
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spurious NN bound states for the cut-off range considered. The resulting NN phase shifts as shown 
in section 3 are significantly improved as compared to the NLO result. We also discuss in this section 
various deuteron properties. It should be mentioned that all these conclusions are based on the type 
of regulator we employ in the Lippmann-Schwinger equation. It cannot be excluded at present that a 
regulator can be constructed that allows for using the large q without leading to deep virtual bound 
states. However, if such regularization exists, it has to look very different than the commonly employed 
regulator functions. 

We then switch to the 3N and 4N systems and briefly demonstrate in section 4 the predictions corre- 
sponding to the NNLO potential. As already stated before, neglecting the 3N forces leads to strong 
deviations from the data. 

The central results of our paper, namely the application of the NNLO* potential to predict 3N and 
4N observables, are presented in section 5. All these results have to be supplemented in the future 
by the inclusion of the three types of topologically different 3N forces which occur at NNLO. This 
additional extensive investigation is left to a forthcoming paper. We summarize briefly in section 6. 

2 Few— Nucleon Forces in Chiral Effective Field Theory 

Starting from the most general chiral invariant effective Hamiltonian density for pions and nucleons 
one can derive nuclear forces by eliminating the pions through a method of unitary transformation [6] . 
Since this transformation acts on the field theoretical Hamiltonian, it leads to an energy-independent 
effective Hamiltonian in the pure miclconic space. The condition for decoupling the purely miclconic 
Fock space states from the ones with pions, a nonlinear decoupling equation, can be linearized by 
introducing a series of orthonormal subspaces with different number of pions leading to an infinite 
set of coupled equations determining the unitary operator. Those equations can be solved recursively. 
Thereby the basic organization principle is a counting scheme in powers of momenta and number of 
pions. We refer to [6] for the detailed steps. Notice also that the relativistic 1/m corrections are 
assumed to be suppressed compared to the 1/A,^ ones, see [4]. Further, we will consider only the 
isospin invariant case in this section. Isospin violating effects can be treated along the lines presented 
in refs.[10],[ll]. The resulting nucleonic potentials are ordered by the power 



where L and Vi are the numbers of external nucleon lines, loops and vertices of type i, respectively. 
Further, the quantity Aj, which defines the dimension of a vertex of type i, is given by 



with di the number of derivatives or M^^ insertions and the number of nucleon lines at the vertex 
i. The inequality Aj > holds true as a consequence of chiral invariance. This leads to > 
for processes with two and more nucleons. One also recognizes that the diagrams with loops are 
suppressed and that (n + l)-micleon forces appear at higher orders than n-nucleon forces. 
Let us now consider first several orders of the NN force. At leading order v = Q (LO) only tree diagrams 
with vertices of Aj = (ttNN vertex with one derivative and two independent four-nucleon contact 
interactions without derivatives) are allowed, see eq. (2.1). Consequently, the LO chiral potential 
is given by the well established one-pion exchange (OPE) and contact forces with the low energy 
constants (LEG) Cs and Ct, as shown in Fig. 1: 




(2.1) 



di + -rii — 2 



(2.2) 



(2.3) 
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Figure 1: Leading order (LO) contributions to the NN potential: one-pion exchange 
and contact diagrams. Graphs which result from the interchange of the two nucleon 
lines are not shown. Solid and dashed lines are nucleons and pions, respectively. 
The heavy dots denote the vertices with Aj = 0. 



.(0) _ ( 9A \ ^ <?1 • g<?2 



^OPEP = -( — ) 'Tl -Ta 



Here p and p' arc the initial and final momenta of the nucleons in the CM frame and q = p' — p. 
Further, M^, gA, and /tt are the pion mass, the axial pion-nucleon coupling constant and the pion 
decay constant, respectively. 

At next-to-leading order (NLO) or ly = 2 there are TPE diagrams with the leading vrNN vertices 
with Aj = according to Fig. 2 and seven contact forces with vertices of Aj = 2 containing two 
derivatives^, see Fig. 3. It should be emphasized at this stage, that the expression (2.1) only allows 
to estimate the order of the corresponding process. It is, however, not possible to read off the precise 
structure of the operators (i.e. the corresponding energy denominators and overall factors) related 
to a particular diagram. This is because the presented figures refer to diagrams within the method 
of unitary transformation and not to ordinary graphs in the old-fashioned perturbation theory. The 
precise operator form of the NLO and NNLO contributions to the 2N and 3N potentials can be found 
in reference [6]. Note also that the graphs 9 and 10 in fig. 2 are not reducible ones in the sense, that 
no energy denominators related to purely nucleonic intermediate states appear in the corresponding 
expressions; see [6] for more details. 

In addition, there are nucleon self-energy contributions and vertex corrections [6] , which renormalize 
the one-pion exchange and contact forces, which we do not show explicitly here. The TPE terms 
shown in Fig. 2 lead to polynomial parts with, in general, infinite coefficients, which renormalize 
various contact interactions, and to finite non-polynomial ones, which are finite and independent of 
the regularization scheme used. The resulting potential reads: 

V^olt = Ci t + ^2 P + (C3 e + Q P)(?i • 02) + iCs \ (<7l + 02) ■ {q X k) 

+ Ce (g ■ ai){q- 02) + C7{k- ai){k ■ 02) , (2.4) 
^^?EP = -^g^^('^){4iV4^(55l-4,i-l) + , 2(23,1 -lO5i-l) + ^0 



*The contact interactions with one insertion of are formally indistinguishable from the four-nucleon operators 
without derivatives and lead to renormalization of the constants Cs, Ct- We will not consider such operators explicitly. 
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Figure 2: First corrections at NLO to the NN potential in the projection formahsm: 
two-pion exchange diagrams. For notations see fig. 1. 



where 

L(,) = iyS|7?toi^5^^±i. (2.5) 

and k = l/2{p' + p). There are seven LEC's Ci to C7 related to contact interactions with two 

derivatives, see fig. 3. 

At that order NLO 3N forces of the topologies shown in Fig. 4 cancel. Note that this cancellation is of 
different type than the one found in time-ordered perturbation theory [12], [13]. To be more precise, 
in that order the contribution of the "irreducible" two-pion (one-pion) exchange diagrams 1-8 (13) 
cancels against the "reducible" two-pion (onc-pion) exchange graphs 9-12 (14, 15). The last graph 
16 in this figure is proportional to the kinetic energy of the nucleons and contributes therefore only at 
higher orders [4]. 

At NNLO (v = 3) there occur new ttttNN vertices with Aj = 1, which contain either two derivatives 
or one insertion and arc parametrized by three constants, denoted in the commonly used notation 
by ci, C3, and C4 (the C2-term does not contribute at this order) [14]. They enter into the TPE NN 
force as shown in Fig. 5 as well as into vertex correction diagrams (not shown), which renormalize 
the OPE, and also into the TPE 3N force shown in Fig. 6. The explicit expression for the two-pion 
exchange NN force at NNLO is^ 

v^?.p - -^i{- ,j:Mi^ + - °3) - + ^)) (2Mi + 



Note that we included here the 1/m corrections, which are formally of the higher order. 
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Figure 3: First corrections to the NN potential: contact diagram at next-to-leading 
order (NLO). The filled diamond denotes seven vertices of Aj = 2 (with two deriva- 
tives). For remaining notations see fig. 1. 



+ 



5127rm/, 



T-^ ((<?! ■ q){d2 ■ q) - q\a, ■ ^2)) i2M^ + q'')A{q) 



327^j4 ("^1 ■ '^2) ■ qW2 ■ q) - q'^i^i ■ 0^2)) 

(p- ' X p) {2Ml + q^)A{q) 
^|^^^(Ti ■ T2) i {a, + a2) ■ (p ' X p) {AMI + q^)A{q) , 



where 



Aiq) = — arctan ^ . (2.6) 

Altogether there arc 9 LEC's at NNLO (and at NLO) related to various contact interactions, which 
have to be fitted by adjusting the NN force to the NN data. The LEC's 01^3^4 which first appear at 
NNLO occur also in ttN scattering and that information should be consistently taken into account. 
The 3N force at NNLO consists of three different topologies as shown in fig. 6. Besides the TPE there 
is a pion exchange between a NN contact force and the third nuclcon and a pure 3N contact force. In 
both cases new vertices of Aj = 1 with unknown constants enter. The precise structure of the chiral 
3NF will be discussed in a forthcoming paper. 

Chiral forces are only valid in a low-momentum region. We enforce this by modifying the above given 
NN force expressions as 

y{v'.v) - /r(p") V^(p",p)/k(p) , (2.7) 
where /i?(p ) is a regulator function. In what follows, we work with the following regulator function: 

rr^v) = eM-pV^') ■ (2.8) 

The power four in the exponent guarantees that the Q^-, and Q^-terms in the potential are 
not affected by the regularization procedure. As already pointed out before, the dependence of the 
low-energy observables on the value of the cut-off A should get weaker with increasing the order v. 
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Figure 4: Leading contributions to the three-nucleon potential at NLO, which can- 
cel: two-pion and one-pion exchange diagrams with the NN contact interaction. 
Graphs which result from the interchange of the nucleon lines and/or from the ap- 
plication of time reversal operation are not shown. In the case of diagram 16, one 
should sum over all possible time orderings. For remaining notations see fig. 1. 



3 Two Nucleons at Next— to— Next— to— Leading Order 

We now turn to the analysis of the 2N system at NNLO. Let us first specify the parameters entering 
the NN potential. The largest uncertainty is related to contact interactions between nucleons. They 



Figure 5: Next-to-next-to-leading order (NNLO) corrections to the NN potential. 
The filled squares denote the vertices with A = 1. For remaining notations see 
figs. 1, 3. 




1 2 3 



Figure 6: Three-nucleon force: TPE, OPE and contact interaction. In the cases 
of diagrams 1 and 2, all possible time orderings should be taken into account. For 
notations see figs. 1 and 5. 



are not restricted by chiral symmetry, but only by the general principles of locality, invariance under 
Lorentz transformations, parity, time-reversal invariance and hermiticity. At NLO and NNLO one 
has to take into account nine independent contact operators contributing to the effective potential: 
two operators without derivatives {vj^l^ in eq. (2.3)) and seven with two derivatives of nucleon fields 
(Krait ™ 6q. (2.4)). The corresponding LECs are fixed by a fit to S- and P-wave phase shifts and to 
ei at low energies. The OPE (F(5?EP in eq. (2.3)) as well as the leading chiral TPE at NLO {V^^p 
in cq. (2.4)) are parameter-free. 

(3) 

As already stressed before, the subleading TPE at NNLO, Vypep ™ ^"l- (^•^)' depends on the LECs 
ci,3,4, which correspond to TTTrNN vertices of dimension A, = 1. Precise numerical values for these 
constants are crucial for various properties of the effective NN interaction as will be discussed below. 
Clearly, the subleading vrvrNN vertices represent an important link between NN scattering and other 
processes, such as vrN scattering. Therefore, ideally, one would like to take their values from the 
analysis of the ttN system, as it was done in [9]. We will now briefly overview the current situation 
concerning the determination of the q's from the vrN system. Several calculations for ttN scattering 
have been performed and published. From the analysis [15] one gets: ci = —0.64 , C3 = —3.90 , C4 = 
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2.25 . Here all values are given in GeV ^. Prom different calculations [15], [16], [17], [18], [19] one 
obtains the following bands for the q's: 

ci = -0.81... - 1.53, C3 = -4.70... -6.19, C4 = 3.25 ... 4.12 . (3.9) 

These bands are also consistent with expectations from resonance saturation, see [16]. Recently, the 
results from a analysis have become available [20]. At this order the S-matrix is sensitive to 
14 LECs (including 01,3,4), which have been fixed from a fit to ttN phase shifts. At this order the 
dimension two LECs acquire a quark mass renormalization. The corresponding shifts are proportional 
to M^. It turns out that different phase shift analyses (PSA) from refs. [21], [22] and [23] lead to 
sizable variations in the actual values of the LECs. A typical fit based on the phases of ref. [22] leads 
to: 

ci = -0.27 ±0.01, C3 = -1.44 ±0.03, C4 = 3.53 ± 0.08 , (3.10) 

where Cj denote the renormalized Cj's. However, using the older Karlsruhe or the VPI phases as input, 
one finds sizable variations in the Cj. Alternatively, one can also keep the Cj at their third order values 
and fit the fourth order corrections separately, see [20]. Due to the uncertainties in the isoscalar 
amplitudes, these constants are not very well determined. The fits could, in principle, be improved 
in the future by including the scattering lengths determined from pionic hydrogen/deuterium. To 
complete the discussion on determination of the q's from the ttN system we would like to stress, that 
numbers consistent with the bands given in eq.(3.9) have been obtained in [24] using IR regularized 
baryon chiral perturbation theory at order and dispersion relations. 

Rentmeester et al. [25] tried to fix the values of the Cj's from an analysis of the pp data, which are of a 
much better quality than the vrN data. In this approach the long-range part of the NN force was taken 
as the sum of the OPE and the chiral TPE (including the NNLO contribution). The NN interaction 
at short distances below some boundary value was parametrized by some artificial energy dependent 
representation. The global fit to the data allowed to pin down the values of the Cj's (and, of course, 
also of the parameters related to the short-range part of the NN force). It turned out that it is not 
possible to fix all three q's in this process because of the strong correlation between these LECs. For 
that reason the constant ci was fixed at the value ci = — 0.76GeV~^ (to obtain a small pion-nucleon 
(T-term of about 40 MeV) and the LECs 03,4 were treated as free parameters. The values of the 03,4: 
C3 = —5.08 GcV^^, C4 = 4.70 GcV~^ determined from the global fit to the pp data are compatible 
with the calculation from the ttN system, see eq. (3.9). Note, however, that this method is not 
directly based on a systematic chiral power counting. 

Having overviewed the current status of the determination of the Cj's from various processes, we are 
now in the position to discuss the corresponding implications for the NN system. First of all, it turns 
out that the numerical values of the Cj's are quite large. Indeed, from a dimensional analysis one 
would expect, for example, the constant C3 to scale like: 

C3~^, (3.11) 

where i is some number of order one. Taking the value C3 = —4.70 from ref. [19] and A^ = Mp = 770 
MeV we end up with i ~ —7.5. Such a large value can be partially explained by the fact that the 
C3,4 are to a large extent saturated by the A-excitation. This implies that a new and smaller scale, 
namely mA — m ~ 293 MeV, enters the values of these constants, see [16]. 

What are the consequences of the large numerical values of the q's for NN scattering? The main 
problem is that the large numerical values of the q's might lead to a slow convergence of the low- 
momentum expansion. To get a feeling of the possible problems one can compare, for instance, the 
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low-momentum matrix elements of, say, the central parts of the TPE at NLO and NNLO. Taking the 
values of the Cj's from the Q^-analysis of the vrN system from ref. [19] 

-4.70 GeV-i , 



ci = -0.81 GeV"^ , C3 = 
as we did in [9] one gets from eqs. (2.4), (2.6): 

. Ml 



C4 = 3.40 GeV~ 



cent, (2) 



TPE 



cent, (3) 



TPE 



q=0 



(4vr/.)2/2 



{l+4g\-8g\) 
6 



~ (ti • r2)(-3.4) GeV" 



(3.12) 



(3.13) 



\2 f2 



(2ci -cs) 10.3 GeV 



-2 



,=0 (4vr/^)^/^ 

Here we neglected all 1/m-corrections. While the order of the matrix element of the potential at 

.cent, (2).,,^,_^)^^M2/(A2/2)with 



\ri-T2} 



NLO agrees with the one expected from dimensional analysis, Vr^-p'E 

£i ~ —0.9, the NNLO matrix element appears to be larger than expected: V^pg' ~ i2M^ / (A^f^) 
where £2 ~ —14.3. Such a deviation from the natural value for £2 of order one does, however, 
not yet necessarily mean a failure of the perturbative expansion, since the potential itself is not 
an observable quantity. To draw a precise conclusion about the convergence properties of the low- 
momentum expansion one should look at the phase shifts, which can be measured directly. Further, up 
to now we only compared the non-polynomial contributions to the potential and omitted all contact 
terms. ^'^ Large numerical values of the low-momentum matrix elements of the Vtpe at NNLO could, 
in principle, be compensated by the corresponding contact terms. However, such a compensation at 
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£^iab [GeV] 
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Figure 7: ^D2 and phase shifts calculated with the CD-Bonn potential. The 
dashed lines show the Born approximation, whereas the solid lines correspond to the 
full solution of the Lippmann-Schwinger equation. The filled triangles are Nijmegen 
PSA results [34] 



NNLO is only possible for S- and P-waves as well as for ei since the contact terms do not contribute 
to D- and higher partial waves at this order. The D- and F-waves may therefore serve as a sensitive 
test of the chiral TPE exchange. The conventional scenario of nuclear forces represented by existing 
OBE models and various phenomenological potentials suggests that the D- and higher partial wave 
NN interactions are weak enough to be treated perturbatively. This is demonstrated in fig. 7 on the 
example of the CD-Bonn potential. Although this observation is confirmed by the smallness of the 
corresponding phase shifts, such a scenario, strictly speaking, does not necessarily need to be realized. 
In fact, the NNLO results can serve as a counter example: with the values of the Cj's from eq. (3.12), 

^^Note that the contact interactions are needed to renormahze the TPE contribution and thus cannot be omitted for 
conceptual reasons. 

"'^"'^This has been suggested by Kaiser et al. in [26], [27]. 
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i^iab [GeV] i?iab [GeV] 

Figure 8: and phase shifts at NNLO using the values of the q's from 

ref. [19]. The dashed hnes show the Born approximation, whereas the sohd hnes 
correspond to the iterated solution with the exponential cut-off A = 1000 MeV. 
The filled triangles are Nijmegen PSA results [34] 



the Born approximation for the S-matrix, for instance, in the D2 partial wave deviates strongly from 
the data already at Siab ^ 100 MeV, see fig. 8. Note that this result is parameter-free and cut-off 
independent.^^ Similar results have been published in ref. [26]. On the other hand, as we showed in 
[9], taking the cut-off of the order of 1 GeV allows for a satisfactory description of all partial waves 
simultaneously. With such a large value of the cut-off, the central TPE potential becomes already 
so strongly attractive that unphysical deeply bound states appear in the D-waves as well as in the 
lower partial waves. Since the potential is very strong (and attractive) and there are no counter terms 
according to the power counting, changing the value of the cut-off clearly leads to strong variation 
of the Z?-wave phase shifts. This is illustrated and discussed in more detail in [9], [28]. Note that 




l^igure 9: and partial waves at NNLO using the values of the Cj's from 
ref. [19]. For notations, see fig. 8. The filled triangles are Nijmegen PSA results [34] 

this problem of the strong cut-off dependence does not show up in lower partial waves, where it is 
compensated by the cut-off dependence of the contact counter terms. In the F-waves, where the 
potential is already sufficiently weak (if the cut-off A is chosen smaller or of the order of 1 GeV) and 
the Born approximation already does a good job, one has no problems with the cut-off dependence 

^■^Since we do not iterate the potential, we do not need to multiply it with the regulating function. Strictly speaking, 
of course, the EFT is only defined with the cut-off procedure which would lead to the results in Born approximation 
being multiplied with an overall factor. For simplicity, we ignore this factor here. 
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as well. This is shown in fig. 9. In spite of this fact one observes sizable deviations for most of the 
F-wave phase shifts from the Nijmegen PSA for energies larger than E'lab ~ 150 MeV [9]. Thus the 

only serious difficulty caused by the large values of the q's in the NNLO analysis of the NN system is 
related to the cut-off dependence of the D-wave phase shifts. 

The large numerical values of the Cj's have also some consequences for three- and more-nucleon 
systems, which will be discussed in detail in the next section. Here we only emphasize that effects 
from the inclusion of the 3N forces are expected to be much larger than in the standard scenario 
of nuclear physics. Note, however, that the separate contributions of the 2N and 3N forces to 3N 
observables cannot be measured experimentally. 

Let us now briefly summarize the consequences of the inclusion of the NNLO TPE with the large 
values of the Cj's taken from the analysis of the ttN system [19]: 

• First of all, including the subleading TPE allows for significant improvement in the description of 

the low-energy observables in the NN system compared to NLO without introducing additional 
parameters, see ref.[9] for more details. The phase shifts are mostly well reproduced. 

• The central part of the potential shows a much stronger attraction than the one found in con- 
ventional models of the NN interaction [26]. As a consequence, one has unphysical deeply bound 
states in the low NN partial waves. 

• The predictions for D-waves depend on the cut-off. The optimal result is obtained for A = 1000 
MeV using the exponential regulator. The potential projected onto the D-waves is strong and 
requires non-perturbative summation via the Lippmann-Schwinger equation. The predictions 
for F-waves deviate from the data at energies larger than Eig^ ^ 150 MeV. In contrast, the 
peripheral waves are well described [26]. 

• One expects large effects from the 3NF. 

Although the NNLO scenario dictated by the large values of the Cj's differs strongly from our expec- 
tations based on the experience with various phenomenological boson-exchange models, one cannot 
exclude this possibility a priori. Indeed, the only serious problem with the large q's is given by the 
strong cut-off dependence of the D-wave phase shifts. However, this will probably not (or only weakly) 
affect chiral predictions for experimentally measured quantities like the cross section, analyzing pow- 
ers, etc. at low energies, where the contribution of the corresponding phases to physical observables is 
rather small. Further, as already discussed in detail in ref.[9], at N^LO it will be cured by dimension 
four contact interactions. Furthermore the failure of the NNLO potential to describe various proper- 
ties in the 3N and 4N systems does not yet indicate a problem, since we have not included the 3NF. 
Because of the calculational difficulties in the treatment of the 3N and 4N systems in the presence 
of deeply bound states it will take some time before all the implications of the chiral EFT at NNLO 
using the large values of the Cj will be explored in detail. These calculations need to be done but will 
require a large amount of computing time. 

Having discussed consequences of the large values of the q's for various properties of few-nucleon 
systems, we can ask ourselves, how confident we are, that the discussed scenario is indeed realized? 
Several comments are in order: 

• First of all, we would like to stress the uncertainty in the determination of the Cj's from vrN 
scattering. The difference between the Cj's from the second and third order analyses of ttN 
scattering is considered to be an effect of third order, i. e. it should be suppressed by one power 
of Q compared to the second order values of the Cj's. For that reason one can equally well 
take the Q^-values of the q's in the NNLO analysis of the NN system, since the Cj's enter 
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only the NNLO and not the NLO contribution to the effective potential. In principle, one can 
also take the values of the Cj from the analysis, which differ from the Cj's by quark mass 

renormalizations of the order M^. Taking different sets of the Cj's from various analyses of the 
ttN system, as described in the beginning of this section, might not cause significant variation 
in description of low-energy observables in the ttN as well as NN systems, but lead to different 
scenarios. 

• It is also possible that including higher order loop effects will reduce the strength of the central 
part of effective NN potential even if the Cj's are numerically large. 

• Finally, already at N^LO one has to include new contact interactions with four derivatives, 
which also contribute in D-waves. These will not only reduce the cut-off dependence of the 
phase shifts, but may also provide additional repulsion and allow to avoid unphysical deeply 
bound states. The work by Entem and Machleidt [29], who constructed a NN potential without 
deeply bound states by a phenomenological extension of the NNLO chiral NN force, may serve 
as an indication of the importance of the higher order contact interactions. To ultimately clarify 
the situation one has to perform a complete analysis of the NN system at order N^LO. 

It is interesting to understand the reason of (possibly) different scenarios in the EFT approach and in 
more phenomenological conventional boson-exchange (BE) models. It has been pointed out in ref. [16] 
that the LECs 03^4 get the dominant contributions from the intermediate A excitation. Also, the a 
and p mesons have been shown to play an important role in the saturation of the q's. In particular, 
the constant ci is completely saturated by the a [16]. Let us now check whether these mechanisms 
of resonance saturation of the Cj's are also realized in the OBE models of the NN interaction. While 
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EFT 
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Figure 10: Exchange of p-meson (wiggly line), which decays into two pions (dashed 
lines) and the corresponding diagrams in EFT (left-hand panel) and OBE models 
(right-hand panel). The shaded blob represents the strong pN form factor in OBE 
models. For remaining notations see figs. 1, 5. 



the resonance saturation of the q's by heavy mesons can, in principle, be interpreted in terms of 
OBE contributions as shown in fig. 10, where the pion loop in the graph in the middle of that figure 
contributes to the form factor of the corresponding heavy meson, the saturation by the A excitation 

^■'To be precise, they included the N^LO contact interactions and allowed for a partial wave dependent cut-off variation. 
Thus, this extension is not an EFT approach. 
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cannot be represented in an appropriate way within the OBE models. Thus, a large portion of the 
subleading chiral TPE is absent in the conventional NN forces. 

A more detailed investigation of the two-pion exchange within the conventional many-boson exchange 
formalism gives rise to a better understanding of the reasons why the intermediate A plays only a 
modest role in the NN interaction. In the Bonn model of ref.[30], which also allows for two-boson 
exchanges, one finds strongly attractive contributions from TPE. Note that this model also takes 
into account A-excitations in the intermediate states. The diagrams with intermediate A-excitations 
have been shown to give the dominant contribution to the uncorrelated TPE. While the TPE model 
successfully describes high angular momentum partial waves, quantitative description of low partial 
waves appears to be impossible. It is even stated in ref. [30] that "the 27r-contribution appears, in 
general, too attractive and a consistent and quantitative description of all phase shifts can never be 
reached" . It was shown that the strongly attractive contribution of the TPE in low partial waves is to 
a large extent canceled by the tt/o diagrams. The authors of ref. [31] came to a similar conclusion. The 




Figure 11: Representation of the vrp exchange diagram within EFT approaches with 
and without explicit A. BEM stays for boson-exchange models. For notations see 
fig. 10. 



more detailed work on correlated vrp exchange has been performed within the conventional formalism 
by Holinde and collaborators, see [32]. In fig. 11 we show one specific example of the irp exchange with 
the corresponding representation in the EFT approach. It is easy to see that the NLO^^ contribution 
to the effective potential from the diagram shown in fig. 11 only leads to renormalization of the 
corresponding LO contact interactions and thus will only influence the S-wave phase shifts. Thus one 
needs to go to higher orders beyond NNLO in the low-momentum expansion to see effects of the Trp 
exchange on the phase shifts in P- and D-waves. The better way to observe the cancellation between 
the vrvr and up exchanges might be to include vector mesons as explicit degrees of freedom in the EFT. 
That would however require a consistent power counting scheme, which has not yet been constructed. 

The study of the TPE within the Bonn model [35] also indicates a very important role of relativistic 
effects for diagrams with intermediate A's. Incorporating relativistic corrections using IR regulated 
covariant baryon CHPT [36] within the EFT formalism has already been shown to reduce the strength 
of the subleading TPE by about 30% [37]. 

Although phenomenological boson-exchange models provide a plausible explanation of the fact that 

^■^Note that if the A-resonance is included explicitly via the "small scale expansion" [33], the strong attractive central 
contribution to the TPE appears already at NLO and not at NNLO. 
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Figure 12: Fits and predictions for the S-waves for nucleon laboratory energies i?iab 
bebw 200 MeV (0.2 GeV). Left/right panels: NLO/NNLO* results. The cut-off is 
chosen between 500 and 600 MeV leading to the band. The filled circles depict the 
Nijmegen PSA results [34]. 



the A~resonance does not play a significant role in NN scattering, additional model independent 
analysis is needed to improve on our understanding of the TPE. In particular, more work on pion- 
nucleon scattering (dispersive versus chiral representation), new dispersive analyses and more precise 
low-energy data are needed to pin down these LECs to the precision required here. 
Motivated by the observed cancellation between the vrvr and vrp exchanges and by the fact that the A 
is not included as an explicit degree of freedom in existing OBE models and is supposed to play only 
a modest role for NN interactions at low energies, we constructed the NNLO* version of the effective 
potential [38] , [39] , in which we basically subtracted the A-contributions from these LECs and allowed 
for some fine tuning. This results in numerically reduced values of the 03^4: 

C3 = -1.15 GeV-^ , C4 = 1.20 GeV"^ . (3.14) 

As a consequence, the attraction of the central potential corresponding to chiral TPE is reduced 
compared to the NNLO calculation of ref. [9]. Differently to the NNLO potential, we also incorporated 
in the NNLO* version the leading isospin violating effect due to the pion mass differences in the OPE. 



We are now in the position to discuss numerical results of the NNLO* potential. First, we make some 
general remarks. For NLO (NNLO*), we fit to the Nijmegen S- and P-wave phases and the ei mixing 
parameter up to £^iab = 50 (100) MeV. These phase shifts at higher energies and for all higher partial 
waves are therefore predictions. Throughout, we show the phase shifts using the exponential regulator 
given in eq.(2.8). We are now able to use the same cut-off range as we did at NLO. Varying the cut-off 
A between 500 and 600 MeV, we find a weakly changing x^/ P^r degree of freedom. Also, for this 
range of the cut-off we do not encounter any unphysical bound state in any partial wave, which is in 
stark contrast to the NNLO results of [9]. We note that one finds an increasing number of such deep 
bound states with increasing cut-off, eventually leading to a limit cycle behavior (for details, see [3]). 
The theoretical predictions at NLO and NNLO* for this cut-off range are indicated as bands in the 
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Figure 13: Fits and predictions for the P-waves and the mixing parameter ei 
for nucleon laboratory energies E'lab befow 200 MeV (0.2 GeV). Left/right pan- 
els: NLO/NNLO* prediction. The cut-off is chosen between 500 and 600 MeV as 
shown by the band. The filled circles depict the Nijmegen PSA results. 

following figures. In most partial waves these bands get thinner when going from NLO to NNLO* and 
are also visibly closer to the data (Nijmegen PSA). This is what one expects from a converging EFT. 
Let us now regard different partial waves. In fig. 12 we show the two S-waves. We find a good 
description at NNLO* up to 200 MeV, which is comparable with (in case of the ^5*0 partial wave 
slightly worse than) the NNLO results shown in ref. [9]. 

Consider next the P-waves and the mixing angle ei shown in fig. 13. The most visible improvement 
from NLO to NNLO* is observed for and ei. We also note that the description of is better 
than in the NNLO case shown in ref. [9]. While the NNLO corrections to the NLO results for the ^Pi, 
3Pi and partial waves (see fig. 5 in ref. [9]) go in the right directions, the observed effects turn out 
to be too large and lead to significant deviations from the data. This is cured in the NNLO* version, 
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Figure 14: Predictions for the D-waves and the mixing parameter 62 for nucleon 
laboratory energies .Ei^b below 200 MeV (0.2 GeV). Left/right panels: NLO/NNLO* 
prediction. The cut-off is chosen between 500 and 600 MeV as shown by the band. 
The filled circles depict the Nijmegen PSA results. 



as can be seen from fig. 13. The NNLO* and NNLO results for the ^Pq partial wave are very similar 
to each other and to the NLO calculation. 

Let us now discuss the D-waves and the mixing angle £2. These are of particular interest since at 
NNLO* no parameters enter and we already discussed the strong cut-off sensitivity found at NNLO. 
As shown in fig. 14, this cut-off sensitivity is sizeably reduced at NNLO* (in comparison to NNLO) 
and one obtains an overall good description of all D-waves up to laboratory energies of about 200 MeV. 
We remark that the important vrvr correlations which are at the heart of the dramatic improvement in 
^1)3 from NLO to NNLO* are still present (as in NNLO) since they are driven by the physics behind 
the LEG ci. Note also the significant improvement for the 62- 

The NNLO* corrections get weaker for F- and higher partial waves. In contrast to the strong NNLO 
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Figure 15: Predictions for the F-waves and the mixing parameter £3 for nucleon 
laboratory energies ^lab below 200 MeV (0.2 GeV). Left/right panels: NLO/NNLO* 
prediction. The cut-off is chosen between 500 and 600 MeV leading to the band. 
The filled circles depict the Nijmegen PSA results. 



effects in the F~waves, which cause significant deviations of the phase shifts for the results of the 
Nijmegen PSA, the NNLO* results can be viewed as small corrections to the NLO calculations, see 
fig. 15. Indeed, in most cases the difference between the NLO and NNLO* predictions is very small. 
The only exception is observed for the ^F/^ partial wave. Here the NNLO* corrections go in the right 
direction but are still not sufficient to reproduce the phase shift appropriately at energies larger than 
50 - 100 MeV. 

The peripheral partial waves (G,H,I, . . .) are mostly well described. Most of these are dominated 
by OPE. However, in very few cases the large values of the Cj's were needed to bring the prediction 
in agreement with the data, see refs.[26, 9]. In the NNLO* potential, the weakened TPE does not 
provide enough strength as e.g. seen in ^Gs, cf. fig. 16. Similar remarks hold for the H and / phase 
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Figure 16: Predictions for the G-waves and the mixing parameter £4 for nucleon 
laboratory energies Ei^y, bebw 200 MeV (0.2 GeV). Left/right panels: NLO/NNLO* 
prediction. The cut-ofF is chosen between 500 and 600 MeV leading to the band. 
The filled circles depict the Nijmegen PSA results. 



shifts; we refrain from showing these here. 

We now turn to the bound state (deuteron) properties. We have not fine-tuned the parameters to 
exactly reproduce the binding energy. It is already described within 2% for the range of cut-offs 
considered here. In table 1 we collect the deuteron properties at NLO and NNLO* (for A = 500 
and 600 MeV) in comparison to the NNLO results (obtained with an exponential regulator with 
A = 1.05 GeV) and the CD-Bonn potential (as one generic high-precision potential). Most deuteron 
properties are well reproduced and improve when going from NLO to NNLO*. We also note that all 
NNLO* predictions (except the one for the quadrupole moment) are between the NLO and NNLO 
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results. The quadrupole moment is only slightly improved at NNLO*, while the NNLO correction 
for this quantity goes in the wrong direction. One, however, still observes a discrepancy of about 7% to 
the experimentally observed value (sec, however, the recent discussion by Phillips [40] why this failure 
is not unexpected). It has also been noted in [41] that fine-tuning the binding energy can slightly 
improve the prediction for Q^,- The NNLO* and NNLO corrections go in the wrong (right) direction 
for the asymptotic D/S ratio r] (the asymptotic S-wave normalization ^5). The improvement for As 
at NNLO* is significant compared to NLO but still leaves space for N'^LO corrections. The same holds 
true for the root-mean-square matter radius r^. We note that the (unobservable) D-state probability 
is reduced as compared to the NNLO result and agrees more with the one found using CD-Bonn 
potential. 





NLO 


NNLO* 


NNLO 


CD-Bonn 


Exp. 


500 McV 


600 McV 


500 MeV 


600 MeV 


Ed [MeV] 


-2.152 


-2.165 


-2.182 


-2.189 


-2.224 


-2.225 


-2.225 


Qd [fm2] 


0.265 


0.266 


0.265 


0.268 


0.262 


0.270 


0.286 




0.0248 


0.0248 


0.0247 


0.0247 


0.0245 


0.0255 


0.0256 


ra [fm] 


1.975 


1.975 


1.970 


1.969 


1.967 


1.966 


1.967 


As [fm-V2] 


0.862 


0.866 


0.871 


0.874 


0.884 


0.885 


0.885 


Pd[%] 


3.17 


3.62 


3.65 


4.52 


6.11 


4.83 





Table 1: Deuteron properties derived from our chiral potential at NLO and NNLO* (for the cut-off 
range considered throughout) compared to the NNLO results of [9], one "realistic" potential and the 
data. Here, is the binding energy, the quadrupole moment, r] the asymptotic D/S ratio, the 
root-mean-square matter radius, ^4^ the strength of the asymptotic S-wave normalization and Pd 
the D-state probability. 

The NNLO* deuteron coordinate space S- and D-wave functions u{r) and w{r), respectively, are 
shown in fig. 17. By construction, they have no nodes and agree quite well with e.g. the CD-Bonn 
wave functions. This lets one expect that the NNLO* potential when applied to the 3N and 4N 
systems gives results closer to calculations based on conventional potentials as does NNLO. We will 
discuss this issue in the following two sections. 

Let us now summarize the presented numerical findings for the 2N low-energy observables. Altogether 
it can be seen that the NNLO* potential leads to results, which are significantly improved compared 
to the NLO ones and allows for a quantitatively rather good description of the np phase shifts up to 
-^lab ~ 200 MeV. While the results for observables at NNLO* and NNLO seem to be of comparable 
quality and in many cases do not significantly differ from each other, these two versions of the chiral 
potential suggest quite different scenarios, as discussed above. It is difficult to give preference to the 
NNLO* or the NNLO version of the chiral potential. In principle, the LECs q should be taken from 
the analysis of ttN scattering and no readjustment should occur, if sufficiently many terms of the 
chiral expansion of the ttN scattering amplitude appear in the TPE potential and the ttN parameters 
are precisely known. However, with the presently available best determinations of these LECs at 
and Q^, one gets a very strong attractive central part of the TPE and, as a consequence, encounters 

^^One should keep in mind, that while the NLO and NNLO* results are given within the theoretical uncertainty, which 
corresponds to a cut-off variation, the results at NNLO are only shown for the optimal choice of the cut-off A = 1050 
MeV. 
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Figure 17: Deuteron wave functions at NLO (left panels) and NNLO* (right panels) 
in coordinate space. The solid (dashed) lines correspond to the cut-off A = 500 
MeV (A = 600 MeV). 

unphysical deeply bound virtual states. Further, they lead to an unconventional balance between 
two- and many-nucleon forces in systems with three (or more) nucleons. On the other hand, boson- 
exchange phenomenology clearly indicates the suppression of contributions with delta intermediate 
states based on cancellations with e.g. tt/j exchanges. Such a scenario is realized in the NNLO* 
potential, which does not lead to unphysical bound states in the NN system for reasonable choices 
of the cut-off. Progress can come from different directions: Further investigations of the vrN system 
at higher orders in chiral expansion as well as new data (eventually combined in dispersion relations) 
may allow for more precise determination of the Cj's, so that one would be able to discriminate the 
physically relevant scenarios of the NN interaction. On the other hand, the final word on the choice of 
regulator is not yet spoken - one may still contemplate the construction of a coordinate-space regulator 
that modifies the TPE at short distances such that no unphysical bound states appear. At present, 
this is only a speculation (we refer to [42] for some related work). Clearly, more work in this direction 
is mandatory. For the time being we consider it legitimate to use the NNLO* potential in applications 
to the 3N and 4N systems. For the sake of completeness, we will, however, briefly discuss in the next 
section NNLO predictions for the 3N system, before we switch to the central issue of this paper and 
present the NNLO* results for 3N and 4N systems. 

4 NNLO Predictions for the 3N System 

As has been shown in [9] the NNLO NN forces describe the Nijmegen NN phase shift values significantly 
better than the NLO ones. We would like to remind the reader that there occur spurious bound states 
in S-, P- and D-waves, as already mentioned in the preceeding section. As a consequence of these 
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deeply bound states, the deuteron wave function at NNLO has nodes below about 1 fm, which are 
not present at NLO (and NNLO*) or using conventional NN forces. In agreement with the correct 
description of the low-energy ^Si Di phase shift parameters those nodes also do not influence the 
low-energy deuteron properties: its binding energy, the asymptotic D/S ratio, the root-mean-square 
matter radius, the asymptotic S-wave normalization constant and the quadrupole moment, which are 
in good to fair agreement with the experimental values. 
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Figure 18: Differential cross section (in [mb/sr]) and analyzing powers Ay, T20 and 
T21 for elastic nd scattering at ii^iab =3 MeV. The solid and dashed lines correspond 
to predictions based on the CD-Bonn and NNLO potentials. The open (filled) 
triangles are pd [43], [44] {nd [45]) data. 



In turning to the 3N system one encounters in the Faddeev formulation NN t-matrices which are 
taken off the energy-shell. The energy argument is 

E2 = E-^q\ (4.1) 

where E is the fixed 3N energy and (3/4m)g^ the varying kinetic energy of the third particle in 
relation to the pair of nucleons interacting via the NN t-matrix. Since q varies between and infinity 
one necessarily hits the spurious bound state energies, which occur as poles of the NN t-operator. 
Physically spoken this has the consequence that the normal 3N bound state is not stable but decays 
into two fragments, a deeply bound spurious NN bound state and a nucleon. In practice this decay 
is rather weak, however, and can be neglected since the physical 3N bound state has little overlap 
with the short ranged spurious NN bound state. In addition, one has to expect that there will be 
spurious 3N bound states at extremely large negative energies in the GeV region. Calculating the 
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3N observables in the presence of deeply bound spurious states in the NN system requires some 
precautions. Of course, the ultimate way to calculate 3N observables in the presence of deeply bound 
NN states would be to treat the poles of the NN t-operator explicitly in the corresponding integral 
equation. The much easier approximate way is to restrict the virtual g-values such that E2 does not 
reach the energies of the spurious bound states, which are in the GeV region. Alternatively, one can 
transform the two-body Hamiltonian in such a manner that the NN phases do not change but the 
spurious bound state energies are moved towards high positive energies, where they cause no technical 
obstacles. This can be achieved for instance by the following simple change of the two-nucleon force 
(starting from the Hamiltonian H = Ho + V) 

V = V + Y,\^iM^i\ (4.2) 

i 

leading to the modified Hamiltonian H = Hq + V with shifted eigenvalues corresponding to spurious 

eigenstates, 

H\^>i) = {Ei + ai)\^i) , (4.3) 

where the aj are sufficiently large positive energies, the spurious bound states and Ei the binding 
energy of the spurious state Note that such a projection does not influence the 2N phase shifts. 

Also the deuteron wave function remains unchanged. 

In this work we do not aim to apply the NNLO potential to the 3N and 4N systems and only want to 
demonstrate that one needs strong 3NFs to describe the data. The approximate methods described 
above are therefore sufficient for our present purpose. We solved the 3N Faddeev equation for the 
triton using V instead of V and this for NNLO. With the cut-off A = 1000 MeV in the NN system 
we found for triton binding energy E = —3.8 MeV. This number turned out to be nearly independent 
of the actual values of Ui (which are of the order of a few GeV). A very close value arises if on sticks 
to the original NN force V at NNLO and restricts the range of q values as mentioned above. One 
has to conclude that this form of the NN force requires strong 3N forces to account for the missing 
binding energy. Notice, however, that these required 3N forces may still be much weaker than the 
corresponding 2N ones. Indeed, we found an expectation value for the potential energy in the triton 
at NNLO of about —172 MeV, which is much larger than the one observed for various high-precision 
potentials of the order —40 to —50 MeV.^^ The corresponding large value for the kinetic energy has, 
in principle, to be expected due to the additional nodes in the deuteron and triton wave functions in 
the short distance range, which are caused by the deeply bound states. 

In view of the results for the triton binding energy one also has to expect that theoretical 3N scattering 
observables based only on the NNLO NN force (i.e. neglecting the 3NFs) will be in conflict with the 
data. This is indeed the case as shown in fig. 18 for a few examples. 

Thus, we conclude that taking into account only the 2N interaction at NNLO and neglecting the 
corresponding 3NFs does not allow for a correct description of the 3N observables. This presumably 

will be corrected by the inclusion of the 3N forces, which because of consistency in the power counting 
has to be taken into account at NNLO. It will be interesting in the future to check this statement 
explicitely. 

5 3N and 4N Predictions with the NNLO* NN Potential 

We use the Faddeev- Yakubovsky scheme to solve for the 3N and 4N bound states and the 3N scattering 
observables as described in [46] , [8] . The calculations are fully converged with respect to the number of 

^^Thus, the missing binding energy of about 4 MeV for the triton to be provided by 3N forces is still much smaller 
compared to the strength of the 2N interaction. 
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partial wave states and standard numerical discretizations. Table 2 shows the results for the 3N and 
4N binding energies using the NLO and the NNLO* NN potentials. Note that for the NLO version the 
numbers slightly different from the ones published in [7] appear since we have now taken into account 
the leading isospin violating effect due to the charged to neutral pion mass difference in the OPE. We 
see a clear reduction of the cut-off dependence in going from NLO to NNLO*, as it is expected from a 
converging EFT. For reasons of comparison, we also display the kinetic energy and the probabilities of 
the various ground state components {S, P, D) in ^He. The resulting binding energies for NNLO* are 



Potential 




^(^He) 


T 


S [%] 


P [%] 


D [%] 


NLO, 500 
NLO, 600 


-8.544 
-7.530 


-29.57 
-23.87 


61.4 
77.6 


94.71 
92.60 


0.07 
0.11 


5.22 
7.29 


NNLO*, 500 
NNLO*, 600 


-8.590 
-8.245 


-29.96 
-27.87 


62.2 
64.9 


93.65 
90.61 


0.10 
0.17 


6.25 
9.22 


AV-18 

CD-Boim 


-7.628 

-8.013 


-24.99 

-27.05 


97.8 
77.2 


85.89 

89.06 


0.35 

0.22 


13.76 

10.72 


exp 


-8.48 


-29.00 











Table 2: Theoretical and ^He binding energies for different cut-offs A at NLO and NNLO* compared 
to the AV-18 and CD-Bonn predictions (point Coulomb interaction pcrturbatively removed), the 
experimental '^H binding energy and the Coulomb corrected '^He binding energy in MeV. The kinetic 
energies T (in MeV) and S, P and D state probabilities for ^He are also shown. 

near the experimental data and larger than the values typically achieved with conventional potentials. 

The results for two representatives, AV18 and CD-Bonn, are also displayed in Table 2. Note that the 
NNLO* results encompass the experimental values, quite in contrast to the realistic potentials. We 
remark, however, that the chiral NN forces employed up to now are for the np system and therefore 
do not yet take all relevant isospin violating effects into account Experience tells us that this leads 
to an unphysical increase in the binding energy of about 200 keV (1 MeV) in (''He). Nevertheless 
in relation to conventional forces one ends up close to the experimental data for and "^Hc using 
the NNLO* NN potential and consequently will need smaller contributions of 3N forces than using 
conventional NN forces. 

For 3N scattering we show in figs. 19-24 clastic Nd scattering observables for laboratory energies of 
3, 10 and 65 MeV, in order, and in fig. 25 Nd break-up cross sections for two arbitrarily selected 
kinematical configurations at -Eiab =13 MeV. In each case the NLO are compared to the NNLO* 
predictions and the ones based on the modern high-precision potentials. Like for the bound state 
energies we find in all cases a much reduced cut-off dependence for NNLO* in comparison to NLO. 
Also, at the highest energy we considered, 65 MeV, one observes now a strong improvement compared 
to the NLO results, which in some cases deviate significantly from the data. We also observe that the 
theoretical uncertainty due to the cut-off variation is sometimes smaller than the spread using the 
various phase equivalent conventional potentials. Note that most of the deviations of the theoretical 
predictions from the pd data in case of the tensor analyzing powers and the differential cross section 
at low energies and at forward angles arc due to the Coulomb pp force [55]. 

^^Such effects can be dealt with in nuclear EFT as discussed e.g. in [11]. 
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Figure 19: Differential cross section for elastic nd scattering in [mb/sr] for ii^iab = 
3, 10,65MeV (top to bottom). Results at NLO (left panel) and NNLO* (middle 
panel). The bands correspond to the range A = 500 to 600 MeV. Results based on 
the high-precision potentials (CD-Bonn, AV-18, Nijm-93, Nijm-1,11) are shown in 
the right panel. Here the bands refer to the spread in predictions using the various 
potentials. The crosses are data: at 3 MeV from [43], at 10 MeV from [47], and 
at 65 MeV from [48]. The circles at 65 MeV are nd data from [49]. 



Let us now take a closer look at the calculated elastic observables. The differential cross section 
at NNLO* agrees well with the data and with the predictions based upon various high-precision 
potentials, cf. fig. 19, and is strongly improved at 65 MeV compared to the NLO results. 
The vector analyzing power of elastic nd scattering at low energies is well known to be underpredicted 
by the standard NN potential models, see fig. 20, right panel, and this remains true even after inclusion 
of the existing 3N forces based on boson exchanges. As reported in ref. [7] and shown in the left panel 
of fig. 20, the NLO predictions at 3 MeV are essentially in agreement with the data, while at 10 MeV 
one even observes a slight overestimation in maximum. The NLO results for Ay at 65 MeV show 
significant deviations from the data. Our predictions at NNLO* are much closer to the results based 
upon the high-precision potentials, i.e. the data are underpredicted at low energies (3 and 10 MeV) 
and reproduced accurately at higher ones (65 MeV), cf. fig. 20. Although some improvement with 
respect to the predictions based upon the high-precision potentials can be seen at 3 and especially at 
10 MeV, the pending puzzle is now back at NNLO*. As pointed out in ref. [56], one possible reason 



25 
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Figure 20: Analyzing power Ay for elastic nd scattering, for £^iab = 3, 10, 65 MeV 
(top to bottom). Data at 3 MeV are from [45] {nd, open circles) and [44] {pd, filled 
circles), at 10 MeV from [50] {nd, open circles) and [47] {pd, filled circles), and at 65 
MeV from [49] {nd, open circles) and [48] {pd, filled circles). For further notations, 
see fig. 19. 



for the significant change of about 20 % in the Ay predictions when going from NLO to NNLO* 
may be the deviations of the np ^Pj phase shifts from the data at NLO. These channels are well 
known to be very important for the nd Ay, see e.g. [8]. In table 3 we demonstrate that these partial 
waves are now much better described at NNLO*. We also remind the reader that in contrast to 
high-precision potential models, which are constructed to perfectly reproduce the NN data below 
the pion production threshold, in EFT one does not aim at a perfect description of the data by 
increasing the phenomenological content of the NN interaction but rather at performing systematic 
order-by-order calculations. At each specific order in the low-energy expansion (in our case chiral 
expansion) one has some theoretical error due to missing higher order terms, which can be estimated. 
Considering our results for Ay at NLO one should therefore keep in mind the level of precision of the 
NLO approximation. Further, since the nd Ay is a very sensitive observable and is strongly affected by 
changing the np ^Pj phase shifts by only few percent, the large uncertainty for this specific observable 
has to be expected in the EFT approach. 

The situation with the deuteron vector analyzing power iTn is very similar to the one with Ay. This 
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Figure 21: Tensor analyzing power iTn for elastic nd scattering, for iiJiab = 
3, 10, 65MeV (top to bottom). The circles are pd data: at 3 MeV from [44], at 
10 MeV from [47], and at 65 MeV from [51]. For further notations, see fig. 19. 



is shown in fig. 21.^^ The NNLO* predictions for the tensor analyzing powers T20 and T21 at 3 and 
10 MeV as well as for T22 at all three energies follow the band made up from the variations among 
the high-precision potentials. Remarkably, our results for T20 and T21 are even significantly closer to 
the data at 65 MeV. 

In case of the specific 3N break-up results shown in fig. 25 the chiral force predictions are equally off 
the data as the predictions of the conventional forces. In case of the upper row the deviations in 
the quasi-free peak to the pd data might be due to Coulomb force effects, whose precise size is still 
unknown. The lower row addresses the space-star anomaly. We underestimate significantly the two 
sets of nd data, which are also far off the pd data. As in the case of elastic scattering observables the 
NNLO* predictions follow the band made up from the various high-precision potentials. Again the 
size of Coulomb force effects is unknown. For more information on these break-up configurations see 
refs. [8], [53]. 

It is also interesting to compare our results to the ones shown in ref. [56], in which the same nd 
scattering observables have been calculated using the phenomenological high-precision extension of 

^^Notice that only pd data exist for this observable. Inclusion of the Coulomb interaction will lead to significant 
underestimation of the iTn [55]. 



27 



Tlab [MeV] 


NLO 
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NPSA 


500 MeV 


600 MeV 
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1 


0.19 


0.19 


0.17 


0.17 


0.18 


0.18 
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1.68 


1.67 


1.59 


1.58 


1.61 


1.63(1) 


10 


3.74 


3.72 


3.62 


3.58 


3.62 


3.65(2) 


25 


8.28 
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8.10 


8.13(5) 
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11.12 


11.01 


10.74 


10.70(9) 


100 


8.27 


8.31 


8.34 


8.43 


8.57 


8.46(11) 


150 


2.52 


2.52 


1.88 


2.15 


3.72 


3.69(14) 


200 


-3.70 


-4.11 


-5.43 


-5.23 


-1.55 


-1.44(17) 


'Pi 














1 


-0.12 


-0.12 


-0.11 


-0.11 


-0.11 


-0.11 


5 


-0.99 


-0.99 


-0.91 


-0.92 


-0.93 


-0.94 


10 


-2.17 


-2.16 


-2.02 


-2.02 


-2.04 


-2.06 


25 


-5.05 


-5.03 


-4.82 


-4.83 


-4.81 


-4.88(1) 


50 


-8.35 


-8.32 


-8.22 


-8.23 


-8.18 


-8.25(2) 


100 


-12.61 


-12.66 


-13.49 


-13.47 


-13.23 


-13.24(3) 


150 


-15.56 


-15.94 


-18.48 


-18.43 


-17.51 


-17.46(5) 
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-17.80 


-18.86 


-23.64 


-23.65 


-21.38 


-21.30(7) 


'P2 














1 


0.02 


0.02 


0.02 


0.02 


0.02 


0.02 
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0.24 


0.24 


0.25 


0.25 


0.26 


0.25 


10 


0.70 


0.70 


0.71 


0.71 


0.72 


0.71 


25 


2.87 


2.89 


2.64 


2.65 


2.60 


2.56(1) 


50 


8.05 


8.29 


6.29 


6.34 


5.93 


5.89(2) 


100 


20.32 


22.60 


11.31 


11.70 


11.01 


10.94(3) 


150 


29.73 


35.97 


12.11 


13.04 


13.98 


13.84(4) 


200 


34.02 


44.30 


9.92 


11.31 


15.66 


15.46(5) 



Table 3: ^Pj np phase shifts at NLO and NNLO* for the smallest and largest values of the cut-off 
compared to the phases based on the CD-Bonn potential [57] and to the Nijmegen PSA [34]. 



the chiral potential by Entem and Machleidt [29]. In fact, our results for these observables show a 
remarkable similarity to the ones presented in this reference, i.e. both predictions agree with the 
calculations based upon the conventional high-precision potential models and with the data in most 
cases and are slightly closer to the data for T20 and T21 at 65 MeV. The only significant differences 
between our results and the ones of ref. [56] are observed for Ay (and iTn) at low energies (3 and 
10 MeV), which are slightly improved in case of the NNLO* version. It is very gratifying to see 
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Figure 22: Tensor analyzing power T20 for elastic nd scattering, for 

3, 10, 65MeV (top to bottom). The circles are pd data: at 3 MeV from [44], 

10 MeV from [47], and at 65 MeV from [51]. For further notations, see fig. 19. 



at 



that at least up to = 65 MeV our NNLO* potential with 11 adjustable parameters works for nd 
scattering equally well as the one of ref.[29] with 46 adjustable parameters. This remarkable agreement 
may serve as a nice demonstration of the power and the advantage of an EFT with consistent power 
counting compared to more phenomenological approaches: performing chiral expansion of the nuclear 
force up to some definite order by inclusion of all relevant diagrams and counter terms allows to 
describe low-energy observables with the same precision regardless of the kind of system the theory 
is applied to (2N, 3N, ...). From the point of view of EFT, it makes not much sense to improve the 
description of the 2N observables alone by a phenomenological extension of the short-range part of 
the NN force. As one can see comparing figs. 19-24 with the corresponding ones of ref. [56], this does 
not lead to an improvement in describing other systems at low energy (i.e. the 3N system). In order to 
reduce the theoretical uncertainty, one should instead go to higher orders, which requires the inclusion 
of 3N, 4N, interactions as well as more pion exchanges in the 2N force. Furthermore, the whole 
concept of developing phenomenological NN potentials, which reproduce the NN data perfectly with 
X^/datum= 1 is in conflict with the general EFT philosophy: at each fixed finite order of the low- 
energy expansion one necessarily has some definite uncertainty in description of observables. Adjusting 
the cut-off parameters in various partial waves to improve the fit to data, as it has been done in ref. [29] , 
is not acceptable from the point of view of pure EFT, where the cut-off dependence of observables 
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Figure 23: Tensor analyzing power T21 for elastic nd scattering, for iiJiab = 
3, 10, 65MeV (top to bottom). The circles are pd data: at 3 MeV from [44], at 
10 MeV from [47], and at 65 MeV from [51]. For further notations, see fig. 19. 



180 



may serve as an estimation of the theoretical error. 

It is now an urgent task to encode the three topologically different 3N forces, which have to be taken 
into account at NNLO (NNLO*) and to determine the corresponding parameters in the 3N system. 
Pioneering studies in [58] indicate that specifically the diagram in fig. [5] of this reference might have 
a chance to solve the Ay puzzle. This extensive work will be dealt with in a forthcoming paper. 



6 Summary 

The concept and the resulting NN forces at LO, NLO and NNLO of xPT have been reviewed. Our ap- 
proach is based on the method of unitary transformation applied to the most general chirally invariant 
Hamiltonian expressed in terms of pion and nucleon fields. This method leads to energy independent 
nuclear forces, a property which is important for the application to more than two nucleon systems. 
The NNLO NN forces driven by the low energy constants 01^3^4 lead to deeply bound unphysical NN 
states in low partial waves if the values 01^3^4 are taken from typical vrN data analysis. While this has 
no negative observable consequences in the NN system, since the spurious NN bound state energies 

^^Note that this statement might not hold true for different regularization schemes. 
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Figure 24: Tensor analyzing power T22 for elastic nd scattering, for iiJiab = 
3, 10, 65MeV (top to bottom). The circles are pd data: at 3 MeV from [44], at 
10 MeV from [47], and at 65 MeV from [51]. For further notations, see fig. 19. 



are outside the realm of validity of xPT, they lead to a scenario for nuclear physics which is quite 
different from the one driven by conventional nuclear forces. First, the central part of the NN poten- 
tial turns out to be much more attractive as it is expected from conventional approaches. Further, 
the predictions for 3N, 4N, binding energies based upon the purely NN forces are much lower, 
far below the experimental values, and 3N scattering observables deviate dramatically from the data. 
Therefore, unlike for conventional NN forces, which to a very large extent describe the data, and 3N 
forces are only needed as a relatively small additional contribution, the 3N force contributions here 
will be very essential. We provided arguments based upon experiences with meson theoretical poten- 
tials supporting the choice of C3.4 constants, which are numerically smaller and where intermediate 
A-contributions are subtracted out. Based on those values we introduced a novel NNLO* NN force 
which describes NN phase shifts with comparable quality as the NNLO one up to about E'lab = 200 
MeV. These NNLO* potential is free of spurious bound states and leads to predictions in the 3N and 
4N systems which are rather close to the ones familiar from conventional high precision NN forces. It 
is now of highest interest to include the 3N forces which should be taken into account at that order 
in xPT. This work is in preparation. 

In contrast to conventional nuclear forces this chiral approach is systematic in the sense of power 
counting and nuclear forces are expected to be constructed in a convergent scheme. Therefore the 
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Figure 25: 3N break-up cross sections in [mb/MeV/sr^] against the arc length S of 
the kinematicahy allowed locus. The pd data (filled circles) [52] show in the upper 
row a peak (in the middle) related to a quasi-free scattering picture. The two sets 
of nd data from [53] (open circles) and [54] (filled squares), upper group, and the pd 
data, lower ones, are at and in the neighborhood of "space-star" configuration. 

step to NNNLO should be performed in order to see whether convergence can be reached and long 
pending problems with conventional forces like the low energy analyzing power Ay can be solved 
without ad hoc assumptions. 
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